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Abstract. We present results from two surveys for radio-loud quasars. 
These are (1) a northern survey reaching Ssghz = 0.2Jy and (2) the 
Parkes QSO survey reaching S2.7GHZ = 0.25Jy, which together cover most 
of the sky. Both surveys use digitised plate material to identify the major- 
ity of the radio sources and CCD images to identify the optically fainter 
sources. The northern survey was targeted specifically at z > 3 QSOs 
and has produced more than 20, one of which, at z = 4.72, is the most 
distant radio source known. This sample is identified to the POSS-I plate 
limit (20 mag) over an area of 3.7sr, or 12 000 sq deg. With the addition 
of CCD identifications, an effective sub area of 2400 sq degrees is 97% 
identified. In the South the Parkes QSO survey now has complete redshift 
information for all 442 flat-spectrum quasars. These cover the redshift 
range z ~ to z ~ 4.5 and show a clear drop-off in space density at 
redshifts above 3. The form of this decline is remarkably similar in form 
to that seen in optically-selected samples of quasars. Since radio emission 
is unaffected by dust this implies that dust has a minimal effect on the 
observed drop-off seen in optical samples. Analysis of our radio-selected 
samples suggests that the decline in space density at high redshift is more 
pronounced for less powerful radio sources. We are planning a new survey 
that will be sensitive to bright quasars with redshifts up to 6. 



1. Introduction 

The mere existence of QSOs and their host galaxies at high redshift places 
interesting constraints on theories of large-scale structure and galaxy forma- 
tion (Haehnelt & Rees 1993, Efstathiou & Rees 1988). However, recent op- 
tical surveys have produced conflicting results for the QSO space density at 
high redshift, due to the complicated selection effects introduced when picking 
out high-redshift QSOs from more numerous galactic stars with similar optical 
colours. -Radio-based selection methods have the advantage of simpler selection 
criteria. Also, since radio emission is unaffected by dust, the intrinsic evolu- 
tion of the QSO population can be observed, unaffected by obscuration. QSOs 
found in radio-based surveys provide an unbiased (by the effects of dust) sam- 




Figure 1. Sky coverage of the three radio-samples that form the basis 
of our surveys. These are the 0.2Jy flat-spectrum sample of Patnaik et 
al (1992) for 5 > 20°, the O.Uy MG-VLA sample (Lawrence et al 1986) 
in the region 0° < 5 < 20° and flat-spectrum sources from the Parkes 
catalogue in the South. Almost the whole sky (excluding the galactic 
plane, shown by the dashed lines) is covered to a depth of 0.25Jy or 
fainter. Only regions that are complete to a well-defined flux limit are 
used for the analysis of quasar evolution. 

pie of damped Ly-a absorption systems that can be used to study chemical 
abundances to the highest observable redshifts. 

2. Two surveys for radio-loud quasars 

Since high-redshift quasars are rare, surveys to reach z > 3 must cover a large 
area. We are carrying out two surveys using three radio samples that cover 
most of the sky to a limit of 0.25Jy or fainter (see Figure [l] and caption). The 
northern survey uses the 0.2Jy flat-spectrum sample of Patnaik et al (1992) 
and the O.Uy MG-VLA sample (Lawrence et al 1986), which together contain 
more than 4400 sources. The southern survey uses the revised Parkes catalogue 
(Wright & Otrupcek, 1990). Since these samples are defined at high frequency 
(5GHz for the northern samples and 2.7GHz for the Parkes sample), they contain 
a high fraction of core-dominated flat-spectrum sources, usually identified with 
quasars. Moreover, spectral index information is available for all the sources, 
allowing us to select flat-spectrum objects. In addition the sources have accurate 
positions from the VLA and/or the Australia telescope, which are necessary for 
making unambiguous optical identifications. 

The first stage of optical identification is carried out using digitised sky- 
survey plates. CCD images are then obtained for the fainter sources. Since the 
northern survey was aimed at finding high-redshift quasars, only the optically 
red, stellar identifications were followed up spectroscopically. For the Parkes 
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Figure 2. Colour-magnitude diagrams for fiat-spectrum sources with 
S > 0.2Jy from the northern survey, (a) APM/POSS-I identifications 
covering an area of 12 000 sq deg. (b) CCD identifications for a 2400 
sq deg sub-area. Stellar sources are shown as filled circles and open 
circles represent galaxies. The z > 3 quasars are shown as squares. 
The two triangles represent red, stellar objects that do not yet have a 
spectrum. The dashed lines show the colour limits used to define the 
spectroscopic sample. 



sample redshifts were obtained for all the stellar identifications. Details of the 
two surveys are given below. 

2.1. The Northern High-redshift Radio-loud QSO Sample. 

POSS-I plates in the E and O bands, scanned by the APM (Automated Plate 
Measuring facility, Cambridge U.K.), were used to identify flat-spectrum (a > 
—0.5, S oc u a ) sources with Ssghz > 0.2Jy to the plate limit of E = 20mag over 
an area of 3.7sr (12 000 sq deg). Their colour magnitude diagram is shown in 
Figure |2|a. Since high-redshift quasars have redder optical colours than their low- 
redshift counterparts due to absorption by intervening Ly-a, only red (O — E > 
1.0), stellar identifications were followed up spectroscopically. This resulted in 
the discovery of 20 z > 3 quasars (Hook 1994, Hook et al 1995, 1996). Another 
5 z > 3 QSOs with radio fluxes O.Uy < Ssghz < 0.2Jy were found in the region 
covered by the MG-VLA sample. 

About 16% of the radio sources remained unidentified after this stage and 
B and R-band CCD identifications are being obtained for these. So far ~ 40 
have been observed and a colour-magnitude diagram for those detected in the 
R-band is shown in Figure ||b. Combining the POSS-I and CCD identifications 
plus observations of these radio sources from the literature (including the identi- 
fication of one source as a QSO with z = 3.82, Vermeulen et al 1996), an effective 



3 



Lya 



NV 



GB1428 



z=4.72 




k 




SilV/OIV] 




o 



6000 



7000 8 
Wavelength /A 



8000 



9000 



Figure 3. The optical spectrum of GB1428, the most distant known 
radio source and the third highest redshift QSO known (Hook & McMa- 
hon 1998). This object was discovered among CCD identifications of 
fiat-spectrum radio sources with S > 0.2Jy. The spectrum shows un- 
usually weak emission lines for a radio- loud quasar, consistent with the 
scenario that it is beamed (Fabian et al 1997). Note that the spectrum 
has not been corrected for the sky absorption feature at ~ 7600 A. 

area of 2400 sq deg is now 97% identified. 10 objects are still not identified to 
a limit of R ~ 23mag. Two new z > 4 QSOs were found during spectroscopy 
of the 10 reddest CCD identifications. One of these, at z = 4.72 is the highest 
redshift radio source known (Hook &; McMahon 1998). Its spectrum is shown 
in Figure |3|. It is also the highest redshift X-ray source known and is probably 
strongly beamed (Fabian et al 1997). 

We are able to compute the luminosity function of radio-loud QSOs at 
high redshift using the z > 3 QSOs discovered among the POSS and CCD 
identifications. The distribution of the QSOs on the radio-luminosity versus 
redshift plane is shown in Figure Q along with flat-spectrum QSOs at z < 3 from 
Dunlop &; Peacock (1990). The binned luminosity function for these two data 
sets is shown in Figure Notice that the luminosity function at z > 3 is lower 
than that at z ~ 2, particularly for sources of lower radio luminosity. 

2.2. The Parkes QSO Sample. 

Initial identifications of flat-spectrum sources (a > —0.4) from the Parkes cata- 
logue were made by Shaver et al (1996a) using COSMOS scans of UKST plates 
to a limit of Bj ~ 22mag. CCD identifications were obtained for the remaining 
optically-fainter sources. Since all the stellar sources have been identified in the 
B band, Shaver et al (1996a) were able to conclude there are no z > 5 QSOs in 
the sample (see the limit in Figure |6|) and hence there must be a drop-off in the 
space density of radio-loud quasars at z > 3. 

The next stage of this study, namely to obtain redshifts for a complete sub- 
sample of 442 stellar identifications (quasars) with S2.7GHZ > 0.25Jy, has recently 
been completed. The most distant QSO found in this sample has z = 4.46 
(Shaver et al 1996b). Figure || shows a preliminary analysis of space density 
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Figure 4. Distribution of z > 3 radio-loud QSOs with S > 0.2Jy from 
the northern sample (solid symbols) on the radio power vs redshift 
plane, plotted with flat-spectrum QSOs with z < 3 from Dunlop & 
Peacock 1990 (open squares). Circles are QSOs from the MG-VLA 
sample; stars represent QSOs from the S > 0.2Jy flat-spectrum sample 
identified on APM/POSS-I scans and triangles are sources from the 
same radio sample identified on CCD images. 



as a function of redshift for objects with radio luminosity greater than Pn m = 
7.2 x 10 26 WHz _1 sr _1 , corresponding to the flux limit of the sample at z = 5 
(for H = SOkms^Mpc- 1 and q = 0.5). The Fi gure shows that the decline in 
space density of luminous Parkes QSOs at z > 3 has a remarkably similar form 
to that seen in optically-selected samples of QSOs. Since the radio sample is 
unaffected by dust this suggests that dust has a minimal effect on the observed 
drop-off in space density of optically-selected QSOs. This also has implications 
for the effect of dust on recent measurements of the star formation rate from 
high-redshift galaxies. We are working on quantifying this result in the context 
of the obscuration models of Fall & Pei (1993). 



3. QSO Evolution & A New Large Survey 

The results from the above surveys suggest that the peak in QSO space density 
occurs at higher redshifts for stronger sources (Fig [?]). Thus a survey to reach 
higher redshifts should concentrate on relatively bright objects and cover a large 
area, rather than reaching very faint limits over a small region. 

We are planning a very large survey for z > 4 QSOs, making use of new 
wide-area datasets. The survey covers 6000 sq deg and is based on flat-spectrum 
sources with S 5G h z > 50m Jy from the 5GHz PMN sample. The 1.4GHz NRAO- 
VLA Sky Survey (NVSS, Condon et al 1994), which now covers the whole sky 
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Figure 5. The luminosity function of flat-spectrum QSOs as deter- 
mined from the data shown in Figure |I[ Notice that the space density 
at z > 3 is lower than at z ~ 2, particularly for weaker sources. Down- 
ward arrows indicate that the error bar extends to zero since there is 
only one point in these bins. 



north of —40°, will provide spectral index information and the accurate positions 
needed to make unambiguous optical identifications. The UKST I-band survey 
of the southern sky is now nearing completion and the plates are being scanned 
at the APM. As in the northern survey described above, candidate high-redshift 
QSOs will be selected based on positional coincidence, red optical colour and 
'stellar' image classification. By using I-band plates this survey is sensitive to 
bright QSOs out to z ~ 6. Extrapolating our previous results to 50mJy, 8 — 24 
radio-loud objects with z > 4 should be found (2 — 6 times the number currently 
known) plus up to 2 with z > 5. 



4. Conclusions 

There are 4 main conclusions from this work. 

• Radio surveys are efficient for finding high-redshift QSOs. There are now 
many tens of z > 3 radio-selected quasars known including GB1428 at 
z = 4.72, the third highest redshift quasar known. 

• The space density of radio-selected QSOs shows a turnover at high-redshift 
(z > 3). Since radio emission is unaffected by dust, this decline at high 
redshift is not an effect of obscuration. 

• The decline in space density seen at high redshift has a similar form to 
that of optically-selected QSOs. 
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Figure 6. Space densities, normalized to z ~ 2 — 3 and plotted as 
a function of redshift, for the Parkes flat-spectrum radio- loud quasars 
with P lim > 7.2 1 26 W Hz" 1 sr" 1 (•). The number of such objects 
found in the individual redshift ranges are 1 (0<z<l), 12 (l<z<2), 
15 (2<z<3), 6 (3<z<4), and 1 (4<z<5); the error bars correspond 
to ±-\/N. The thick curve is a cubic fit to these data. The upper 
limit shown in the redshift range 5 < z < 7 is taken from Shaver et al 
(1996a). For comparison, similarly normalized space densities are also 
shown for the optically-selected quasar samples of Warren et al. (1994) 
(□), Schmidt et al. (1995) (o), and Hawkins & Veron (1996) (.)• The 
thin lines represent luminosity functions from Boyle (1991) and Hewett 
et al. (1993) used by Warren et al. and Schmidt et al. respectively as 
lower-redshift continuations of their high-redshift space densities. 

• The two previous results lead to the conclusion that the effect of dust on 
the observed QSO drop-off is minimal. 

References 

Boyle B.J. in Relativistic Astrophysics, Cosmology and Fundamental Physics 
(eds Barrow J.D. et al), Ann. NY Acad Sci. 1991, 647, 14 

Condon, J. J. et al., 1994 ADASS III, ASP, ed. Crabtree et al., 155 

Dunlop J.S., Peacock J.A., 1990, MNRAS, 247, 19 

Fabian A.C., Brandt W.N., McMahon R.G., Hook I.M., 1997, MNRAS, 291, L5 
Fall S.M., Pei Y., 1993, ApJ, 402, 479 
Efstathiou G., Rees, M. 1988, MNRAS, 230, 5 
Haehnelt M., Rees M. 1993 MNRAS, 263, 168 
Hawkins M.R.S, Veron P., 1996, MNRAS, 281, 348 

7 




CO 
I 




2 1 2 3 4 5 6 7 



Figure 7. Space density as a function of redshift for sources with 
various radio luminosity P. The top three curves are from the northern 
sample combined with data from Dunlop & Peacock (1990) -see caption 
of Figure ||. The lower curve is from the Parkes sample (see Fig. ||) 
and is plotted with arbitary normalisation. The dotted lines show the 
extent of the redshift bin used to compute each point. 
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